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Abstract: We study holographic models of the QCD temperature-chemical potential 
phase diagram based on the D3/D7 system with chiral symmetry breaking. The baryonic 
phase may be included through linked D5-D7 systems. In a previous analysis of a model 
with a running gauge coupling a baryonic phase was shown to exist to arbitrarily large 
chemical potential. Here we explore this phase in a more generic phenomenological setting 
with a step function dilaton profile. The change in dilaton generates a linear confining qq 
potential and opposes the screening effect of temperature. We show that the persistence 
of the baryonic phase depends on the step size and that QCD-like phase diagrams can 
be described. The baryonic phase's existence is qualitatively linked to the existence of 
confinement in Wilson loop computations in the background. 
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1 Introduction 

Holography [1, 2] has provided a new and powerful description of strongly coupled gauge 
theories. Although a true holographic description of QCD has not been found, many 
insights have nevertheless been obtained through understanding confinement [3, 4], chiral 
symmetry breaking [5-7], phenomenological models of mesons such as AdS/QCD [8, 9], 
and the transport properties of non-abelian plasmas [10]. For further examples, see [11, 12] 
and references therein. Recently [13, 14] we have been interested in finding a holographic 
description of the QCD phase diagram and we will continue those studies here. 

The QCD phase diagram is characterized by the quark and baryon densities and the 
chiral condensate. It is therefore sensible to begin with a simple holographic model that 
encodes the physics of quarks and chiral symmetry breaking. The D3/D7 system is the 
simplest such example [11, 15-18]. It describes the SU(iV c ) M = 4 super Yang-Mills theory 
with Nf quark hypermultiplets. In the quenched approximation the theory is conformal 
and on the gravity side is described by probe D7 branes in AdS§ x 5 5 . The theory is 3+1 
dimensional at all energy scales and has a conformal UV in which the identification of the 
operator matching between the field theory and the gravity description is clean. 

The simplest example of chiral symmetry breaking in the D3/D7 system is found by 
imposing a background magnetic field[19]. In [20] some of the current authors explored 
the temperature-chemical potential phase diagram of this top down system which was sur- 
prisingly complex. The phase diagram is shown in Fig la The chiral restoration transition 
was found to be first order with temperature and second order with density. A critical 
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(a) T-[i phase diagram of D3/D7 system with a 
background magnetic field [20] . The inset diagrams 
are representative probe brane embeddings (dotted 
lines), where a black disk represents a black hole. 
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(c) The phase diagram of the massless ax- 
ion/dilaton gauge theory in [14] showing a baryonic 
phase. 
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(d) The QCD-like phase diagram for a step func- 
tion dilaton computed in this paper: A = 10, A = 



1.715 



Figure 1. Example holographic phase diagrams (see the appropriate references for conventions): 
red lines are first order transitions and blue second order. Phases are labelled by whether chiral 
symmetry is broken or not and whether there is a quark density, d, present. In each case the vertical 
axis is the size of the black hole in the geometry which measures the temperature T. 



point lies between these regimes. In addition there is an extra transition associated with 
the formation of density and the mesons of the system melting into the background plasma 
[21, 22]. In places in phase space this merges with the chiral transition but in other places 
it separates and can be either first or second order. Such a phase with a quark density but 
chiral symmetry breaking could potentially exist in QCD. On the gravity side the phases 
are the associated pictograms shown in Fig 1: the D7 can lie flat, fall on to the black hole 
which provides temperature, or bend in the space to avoid the origin. A number of other 
top-down analyses of related phase diagrams can be found in [14, 23-30]. 

Although the phase diagram of the magnetic field case is interesting it does not over- 
lap with expectations in QCD where the finite temperature transition is believed to be 
second order and the finite density transition first order. In [13] we introduced a more 
phenomenological model in which we treated the running of the gauge coupling/dilaton 
profile as a free function. This was inspired by the fact that the DBI action for the probe 
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with the magnetic field present is equivalent to the same theory with a particular choice 
of unbackreacted dilaton. Using an ansatz which steps between a conformal UV and a 
conformal IR with a larger coupling value we could reproduce phase diagrams like that for 
the magnetic field case but also a wider set of phase diagrams [13]. To describe one that 
looked like the QCD phase diagram we used a further phenomenological freedom. The 
D3/D7 system if backreacted would display an 50(4) x50(2) symmetry in the directions 
transverse to the D3 rather than the SO (6) symmetry of AdSs x 5 5 space. We therefore 
allowed ourselves to introduce a parameter breaking the symmetry in that fashion in the 
emblackening factor of the black hole providing temperature. With this extra feedom we 
could produce the phase diagram in Fig lb which maps more closely to that expected in 
QCD. The high T and low // transition can therefore be fitted well. Here we turn to the 
phase structure in the low T, high fi regime. 

In [14] the phase diagram was studied for the holographic description of a theory with 
a running dilaton. The background has an induced vev for both TrF 2 and TrFF; in the 
gravity dual these scalar solutions satisfy the equations of motion whilst leaving a pure 
AdS background. Quarks are again introduced using a D7 probe. The phase diagram is 
shown in Fig lc and shows the same three phases as the cases already discussed. The 
chiral transition is first order throughout the plane. The extra key component of the 
analysis in [14] was to note that in the geometry with a running dilaton a baryonic phase 
was also present. A baryon vertex is described by a D5 brane wrapped on the 5 5 of the 
AdSs x 5 5 space [31-36]. In the pure M = 4 theory such vertices shrink to zero size. 
However in the running dilaton geometry the large IR value of the dilaton stabilizes the D5 
embedding. Solutions exist that link the D5 to the D7 brane embedding with a balancing 
force condition. These configurations describe the gauge theory with finite baryon density 
rather than finite quark density. The D7 bends off axis to meet the D5 and so the phase 
has chiral symmetry broken. This phase sets in at a particular finite chemical potential and 
then persists to infinite chemical potential (as shown in Fig lc) which is certainly unlike 
the equivalent phase in QCD (the phase diagram is very similar to that found in the Sakai 
Sugimoto model in [37]) . 

The goal of this paper is to combine the approaches of these papers discussed. We will 
investigate the baryonic phase in the model with phenomenological freedom in the running 
coupling. As we had hoped, we will show that these models can possess a baryonic phase 
and that for regions of parameter space it can be in the finite region one would expect it to 
be in for QCD. Our dilaton ansatz transitions between two conformal regions. In each of 
those conformal regimes, as in the N = 2 Karch Katz theory, the baryon vertices are free 
to shrink. Around the transition radius though there is an extra cost to the D5s shrinking 
further since they must encounter the larger dilaton value within. The result is a set of 
stable D5 configurations when the change in dilaton between the two conformal regimes is 
sufficiently large. These configurations are present for low temperatures where the black 
hole horizon has not eaten the D5s but not present at large temperatures. The D7 branes 
that introduce the quark fields can end on the D5 baryon vertex when they are available. 
The D7 branes also prefer not to enter the interior of the geometry with a large dilaton 
value. Thus if the dilaton is sufficiently large inside the D5 radius they will prefer to join 
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Table 1. The brane profile showing the coordinates we use: the background D3, the compact D5 
and the probe D7 

the D5. For very large interior dilaton values this occurs for all values of chemical potential 
and the baryonic phase exists out to infinite chemical potential as in the analysis of [14] (in 
that model the dilaton actually diverges in the IR). However we find that if the change in 
the dilaton is more modest across the transition then whether the D7 energetically prefer 
to join the D5 or end on the black hole horizon must be calculated and depends on the 
parameters. In particular we do find ranges of parameter space where the baryonic phase 
lives in the regime shown in Fig Id corresponding to the region expected in QCD. This is 
the main result of the paper which we will explore in detail. 

Finally we briefly explore the link between confinement and the baryonic phase. The 
quark anti-quark potential can be computed by drooping a string into the AdS interior 
[3, 4]. The step dilaton profile induces a linear potential between a quark and an anti- 
quark when the step size is sufficiently large. The transition from a Coulomb potential to 
linear confinement qualitatively happens at the same value of the step size that induces a 
baryonic phase in the theory at finite chemical potential as one might expect. 

2 The Holographic Model 

We seek to describe a strongly coupled gauge theory with quenched quark fields using the 
D3/probe-D7 system. A running coupling will be imposed on the theory through a radial 
dilaton profile. Finite temperature will be included through an AdS-Schwarzchild geometry 
and baryon number chemical potential via a gauge field in the probe D7 DBI action. The 
final ingredient will be to also allow D5 configurations wrapped on the S 5 of the geometry 
to represent baryon density. We review each of these steps in turn: 

2.1 The Background 

Our geometry will be the AdS-Schwarzschild black hole metric in Einstein frame which 
reads 

ds 2 = (-fdt 2 + d*1) + ^ + R 2 dn 2 , (2.1) 
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where R is the AdS radius {R = ^ng s N c a ) and / is the emblackening factor 



i 



f(r) = l-JL, r H = 7TR 2 T, (2.2) 



with the temperature T. There is also the self-dual RR five-form 

G {5) = (l + *)dt AdfAd(fa). (2.3) 

We will impose an unbackreacted dilaton profile in the AdS-Schwarzschild geometry to 
allow us to explore a range of gauge coupling runnings in the gauge theory. In particular 
we will use the simple step-function form 

= A + 1 - A tanh(r - A) , (2.4) 

where A is the amount of e^'s change at the position A. 
The rescaling of the radius (r) coordinate [5] 



dr dw 1 / „ r~. T Iw^ + w 4 



with \/2wh = thi makes the M 6 structure in the transverse space explicit: 

dfi2 = ^(~9tdt 2 + g x dx 2 ) + ^(dw 2 + dJ2|) , (2.6) 

where 



{w A — wfj) 2 w 4 + Wjj 

w 4 (w 4 + WJj) 



9t — ZI7Z4 , I 4 \ ' 5x — —4 • (^-'J 



2.2 D7 flavor brane and quark phase 

We include a small number (Nf) of quark flavours by placing probe D7 branes into the 
background we have constructed. For this purpose, it is convenient to parameterize M 6 so 
that M 4 x IR 2 is explicit (see the middle of Table 1): 

dw 2 + w 2 dQ 2 5 = dp 2 + p 2 dfll + dL 2 + L 2 d0 2 , (2.8) 

where the D7 brane lies in the (t, x, p, f^) directions and L and p are related to w, 9 



w 



\J p 2 + L 2 , p = wcos9, L = w sin 0, (2-9) 



where < 9 < ir/2 and is different from < 9 < ir in the bottom of Table 1. 

The action for a D7 brane is given by the DBI action without the Wess-Zumino term 

S D7 = -T 7 J d 8 Ze^-det(P[GU + (W)^) (2-10) 
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where T D7 = g ^ (2 ^ )7a , 4 and 



A + 1 - A tanh 



\p 2 + L 2 ) 2 + w 4 H 
p 2 + L 2 



(2.11) 



We add a chemical potential to our system by allowing At = 2irAt(p) ^ 0, and assume that 
the D7 brane embedding profile is none trivial on p only, L = L(p). The action becomes: 



Sd7 = Tm / d xdpC D7 



C D7 = -e*p\lg t g3(l + L> 2 )- 9 4-A? 



(2.12) 



where Tdt = T^NfU^. 

Since At is a cyclic coordinate, there is the conserved charge defined by 
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in terms of which, the Wick-rotated (Euclidean, t 
(A't ->■ d) is defined by 



(2.13) 

-He ) Legendre transformed action 



E LT — / i E 1 

S D7 = T D7 / dt E dx / dpC D ' 7 
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(2.14) 



+ wjj ef 



where the following relation obtained in (2.13) is used. 



A' = d W 
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H 



ur + w 



1 + L' 2 



w A +w 4 H 



p b + 



(2.15) 



w*+wt r 



We solve the equation of motion for L(p) minimizing the action (2.14) for a given d. 
The equation is second order and we impose two boundary conditions: L(oo) = m and 
the IR (p — )• pn) regularity. Asymptotically at large p the embedding field L takes the 
form L = m + cj p 2 + ... with m and c proportional to the quark mass and condensate 
bilinear. By plugging the classical solution L(p) into (2.15), we also get the solution A t , for 
which the large p behaviour is At = p — jt- + ••• with p and d proportional to the chemical 
potential and quark density. To compute the chemical potential, p, we simply integrate 
the right hand side of (2.15), i.e. 



~ w 
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Pn 



dpd 
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1 + L' 2 



n O J n> r 

^ ^ w 4 +w 4 H e 



d 2 



(2.16) 
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where At(pn) = by a regularity condition at horizon. 

The field theory free energy density (J-) is holographically identified as 




oo 



T = T D1 



dpC 



E,LT 
D7 



on— shell 



(2.17) 



where C D ' 7 is written in (2.14). The field theory grand potential density (Q) is identified 
with the Euclidean original action before the Legendre transformation: 




PH 




In both cases, we added the counterterm 



i 



to renormalize the actions. 



l 



cut — off 



2.3 D5 Baryon Vertex 

The D7 world volume gauge field (dual to the chemical potential) has to be sourced by 
strings (quarks in the gauge theory). The string endpoints on the D7 world volume are 
point charges and the world volume gauge field At couples to this point source. The other 
end of the strings must end on either the black hole horizon or another brane. In the 
previous section, we considered the first case: the D7 brane touches the black hole horizon 
and the source for the gauge field are behind the horizon. Now we turn to the latter case. 

Baryons, bound states of N c quarks, are described in AdSs x S 5 by a baryon vertex, a 
D5 brane wrapped on the five sphere with N c fundamental strings attached to it. Therefore, 
it is natural to consider a baryon vertex as a source of strings coupling to the At world 
volume field of the D7 brane. 

For a D5 brane baryon vertex configuration, it is more useful to parameterize M 6 as 



rather than (2.8), since the D5 brane will lie in the (t, O4) directions with a non-trivial 
profile £(0). See the bottom of Table 1. We renamed w to £, to make clear that the 



the D7 brane case. Furthermore, here is also different from 9 in the D7 brane case as 
noted below (2.9). The plots in the following section show the D7 and D5 brane embeddings 
simultaneously in one plot, which is not, strictly speaking, correct because of these different 
coordinate systems. The superposition of the two pictures, which only match where the 
two branes join, is though helpful to understand the solutions. 

1 Since we consider only the massless case, L(oo) — > 0, the counter term is simple. In general, the counter 
term is a function of m and c. 



dw 2 + w 2 dn 2 5 = d£ 2 + f (d9 2 + sin0 2 dfi2) , 



(2.19) 



radial coordinate is a function of 6 here, which is different from w(p) = y L(p) 2 + p 1 in 
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The action for the D5 baryon vertex in the string frame is given by the DBI and the 
Wess-Zumino term: 

S D5 = -T D5 J d 6 £e-^-det(P[GU + (2vra')F ab ) + T m J d 6 U(i) A G {5) (2.20) 

where Tqs = p — ^ 5n , 3 and A^ is the world volume gauge field one-form on the D5 
brane, which is different from At introduced on the D7 brane. In (2.19) the dilaton reads 



A + 1 - A tanh 



'e 4 + < _ A 



(2.21) 



Note that the Wess-Zumino term contributes to the action as well due to the coupling of 
the worldvolume gauge field A(\) with the background five- form Gr$\ (2.3). A nontrivial 
temporal gauge field At(9) couples to N c charge on the D5 brane. With assumptions 
£ = £(0) and At = At(0), the action is given by: 



Sd5 = T D5 / dtd9C D 



where Td§ = R^Q^Tds- Note that we extract i? 4 to make Cd5 dimensionless. 
The equation of motion for the gauge field At(0) takes the form: 

d e D{9) = 4sin# 2 , (2.22) 

where the conjugate momentum D of At is defined by 

D{9) = ^ = A ' Sin4g (2.23) 



The general solution of (2.22) is: 

D(9) = -(vk-9+ -sin0cos6> + sin0 3 cos0), (2.24) 

where < v < 1 is the integration constant and it is related to the number of fundamental 
strings attached to each pole [32]. For our purposes we choose v = which is the case 
that all the fundamental strings emerge from only one pole of the baryon vertex which we 
choose to be 9 = ir. Using (2.23) we can rewrite our Euclidean baryon vertex Lagrangian 
as: 



£%s = W g4^4 + ^ + ? 2 ) ^D{9Y + sin 8 9 , (2.25) 
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where D(9) is defined in (2.24) with v = 0. This is for one Bayron vertex. The holographic 
energy density of many non- interacting free baryon vertex system at finite density n B is 



where hb = and n q is the quark density. 



(2.26) 



2.4 Baryon phase: D7 + D5 branes 

The baryon phase is constructed by connecting a D7 flavour brane and D5 baryon vertices 
by strings between them. It can be shown that the strings' tension is so strong that they 
tend to shrink to a point [14], which makes the D7 flavor brane and D5 baryon vertices 
meet at a point. Therefore, ignoring strings, we start with the D7 - ne D5 combined 
system. Its free energy density is 



•F B = T07 



, r E,LT . n q7f^ 

iv c JO 
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= Tr > i(J : D7 + Fds) 



(2.27) 

(2.28) 
(2.29) 



where 



d since n q is identified with ^ sji t {oo) ^ 3 * s ^ e ^ ree dimensional volume). 



2na'T D7 

Let us consider, at a fixed finite density d, the configuration of the D7 brane with a fixed 
boundary value L(oo) = m and the D5 brane with a fixed £(0) = £o- The two brane 
embeddings have to meet at p = and 9 = tt, i.e. L(0) = ^(tt) = wo. There are infinitely 
many configurations satisfying this condition, parameterized by wo. To find out the lowest 
energy configuration, we vary the total free energy. 



5F B 



dp(EOM L )5L + m 
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(2.30) 
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where the EOM^ and EOM^ are the equations of motion of L and £ respectively. They 
vanish since we consider only the solutions of the equation piecewise. <5L(oo) = <5£(0) = 
by our boundary condition. At the matching point, 5L = 5£, so the condition is reduced 
to 



dC 



E.LT 
D7 



dU 



p=0 



2 fC% 5 
3vr d? 



L'(0) 



£(tt) 



(2.31) 



which is called a force balancing condition [14, 36]. 
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Figure 2. Examples of the baryon vertex solutions and their behaviour as we increase density - 
for these parameters there are two vertices for each value of density. The energetically prefered one 
is always that with the largest radius. The red lines corresponds to d = 0.01, blue d = 5 and green 
d = 1000. Parameters w H = 0.1, A = 10, A = 1.715. 

In the baryon phase, the chemical potential has an extra contribution from D5 branes. 

A = ^~^f = fa + f f dec*, . (2.32) 

where poi means (2.16) now integrated from the D5/D7 join to infinity. Note the extra 
term from the D5 originates from (2.28) where there is an extra coefficient of d. 

However, both in the quark phase and the baryon phase, the grand potential is com- 
puted as a D7 brane Euclidean on-shell action, which is written explicitly in (2.18). The 
D5 brane action does not explicitly contribute in the baryon phase because 

&B = — — = J~b — pd 



T m + r m - U m + — &ec E m \ d 

Fdi — fi-Did . 



However, the D5 brane still implicitly contributes by changing the classical embedding 
solution. In the following section, we will drop the subscripts B and D7. 



3 Phase diagram in grand canonical ensemble 

We will explore the phase diagram of the model with massless quarks. The phase diagram, 
neglecting the baryon vertex, was explored in [13]. The scale A is the only conformal 
symmetry breaking scale in the model and so its value can be scaled (for numerical work 
we take A = 1.715 to match previous work). The phase structure depends on the value of 
the parameter A that determines how much the dilaton changes between the UV and the 
IR. 

For each point in the phase diagram one numerically seeks all possible D7 embedding 
solutions that asymptote to L = at large p. There is always the flat solution L(p) = 0. 
There can also be "Minkowski" solutions that end at p = with L'(0) = 0. Finally there 
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(a) w H = 0.1 (b) w H = 0.345 




(c) w H = 0.8 (d) w H = 1.5 



Figure 3. Behavour of the baryon vertex with temperature for fixed density. There is always 
one baryon vertex for the temperatures considered and typically two as shown. The energetically 
prefered one is always the one with the largest radius. Parameters d = 5, A = 10, A = 1.715. 

can be embeddings that end on the black hole horizon. It is convenient to fix the density d 
and then determine (i from (2.16). The grand potentials of the solutions are then compared 
to determine the prefered phase. 

For A < 2.1 the model is always in the chirally symmetric phase. The massless D7 
embedding is L(p) = at all temperatures and density. To induce chiral symmetry breaking 
the cost of entering the interior volume with the enlarged dilaton must outweigh the cost of 
bending off axis. For A > 2.1 such chiral symmetry breaking is preferred at low T and \i. 
In the parameter regime 2.1 < A < 15 the phase diagram divides into a chiral symmetric 
phase at high T, \x and a region with chiral symmetry breaking at low T, [i. The transition 
between is first order throughout the phase diagram. For A > 15 D7 embeddings that 
end on the black hole play a role in the phase structure at high \x and the phase diagram 
mutates to the form shown in Fig la [13]. The extra phase has non-zero quark number 
and chiral symmetry breaking. The transitions between the phases at high density become 
second order (the transition with T at \x = remains first order no matter how large A). 

In addition to these embeddings we can now seek linked D5-D7 configurations as well 
with the same massless asymptotic boundary condition on the D7. Some examples of these 
configurations are shown in Fig 2 and Fig 3. The D5s are all stabilized for radii of order A. 

Generally the number of baryon vertex solutions (linked D5/D7 solutions with massless 
boundary conditions on the D7) vary with the parameters A, d and T. For some parameter 
sets we don't have any baryon vertex solutions (for small A's and densities). In Fig 2 we 
show an example of a parameter set where we have two baryon vertex solutions. Regardless 
of the number of baryon vertex solutions, we have always found that the energetically 



- 11 - 




(e) w H = 0.1 (f) w H = 0.8 



Figure 4. Plots showing the variation of a number of order parameters through the transition 
regions. The colour coding is - Grccmflat embedding, Black: Minkowski embedding, Blue: large 
radius baryon vertex, Orange: small radius baryon vertex, Red: black hole embedding. Parameters: 
A = 10, A = 1.715, r = 1. Top: the grand canonical potential vs chemical potential, Middle: d-fi 
plot, Bottom: condensate c vs /x. 

favoured solution is the one with the largest radius at the south pole. As we increase A 
the number of solutions increases and the D5 embeddings grow in radius. The baryon 
vertex solutions also change with density. Starting from very small density and increasing 
density the baryon vertex increases its radius rapidly in the beginning before asymptoting 
to a slower growth for larger densities. At fixed density, increasing temperature typically 
generates a smaller baryon vertex due to the black hole attraction (Fig 3). 

Our job now is to again compare the Grand Potential energies of all these configura- 
tions. Some examples of that process are shown in Fig 4(a,b). The points where transitions 
occur can be read from the energies and then checked against transitions in order parame- 
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5 10 15 



(c) A = 20 

Figure 5. Phase diagrams as a function of the parameter A in the dilaton profile. The temperature 
T is in units of wh- Transitions marked in blue are first order, those in red second order. The 
phases are labelled by whether chiral symmetry is broken and whether there is a baryon or quark 
density. 

ters of the model. For example in Fig 4(c,d) we plot the quark density against the chemical 
potential where the transitions and their orders can also be seen. The parameter c deter- 
mining the quark condensate is also shown in Fig 4(e,f). 

After this work one can construct the full phase diagram which is shown in Fig 5 for 
three choices of the parameter A. For small A (but still above the threshold for there to 
be chiral symmetry breaking) the baryonic phase plays no role. The cost of entering the 
high A area is not large enough to discourage the flat embedding. 

At intermediate A we see the baryonic phase enters at intermediate chemical potential 
and low T. Fig 5b looks similar to expectations in QCD for the position of the baryonic 
phase and this is the most significant result we present. Note that the transition to the 
baryonic phase from the vacuum phase is second order. In QCD it is expected to be first 
order at low temperature due to the interactions between the nucleons. At low \x the 
phase in QCD can be thought of as sparse liquid droplets of nucleons (ie nuclei) in the 
vacuum. When [i is sufficient to fill space with nucleons the droplets rearrange themselves 
into the hadron gas. The entropy change associated with the internuclear interactions in 
this rearrangement generates the first order behaviour. In our analysis such interactions 
are neglected so not surprisingly a simple second order transition results. Recent work on 
computing inter-baryon forces using holography are summarized in [38]. 

Finally for large A the baryonic phase is stable out to infinite chemical potential - 
the cost of entering the large A region is so great that the flat embedding is always less 
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Figure 6. The canonical ensemble phase diagram for the case A = 10. The temperature T is in 
units of wh- Transitions marked in blue are hrst order, those in red second order. The phases are 
labelled by whether chiral symmetry is broken and whether there is a baryon or quark density. 

prefered than a D5 ending embedding. This matches the model of [14] in which the dilaton 
diverges at small radius (and also the Sakai Sugimoto phase diagram in [37]). The A = 20 
diagram also displays a small phase region in which black hole embeddings are dominant 
corresponding to a chiral symmetry broken phase but with quark density. 

4 The canonical ensemble 

The phase diagram of our theory shows some additional structure in the canonical ensemble 
that is worth mentioning. We will work just in the case A = 10. Here we fix the density 
d rather than the chemical potential /i and we minimize the free energy rather than the 
grand potential. We display the phase diagram in Fig 6. The chiral symmetry breaking 
phase with d = in the grand canonical ensemble lies entirely on the T axis in this plot. 
Interestingly though the canonical ensemble phase diagram includes a region with black 
hole embeddings (the chiral symmetry breaking, d ^ phase) which is absent in the grand 
canonical ensemble. This is not a discrepancy. At the first order transitions in the grand 
canonical ensemble there is a jump in all order parameters including the quark density d 
- there must therefore be regions of the canonical ensemble phase diagram that are not 
present in the grand canonical ensemble and indeed are energetically unstable. On any 
given fixed T line the first order transition from either the Minkowski embedding or the 
baryon phase to the symmetric phase "leaps" over the black hole embedding phase of the 
canonical ensemble. As we saw in Fig 4(c) for A = 20 this third phase does play a small 
role in the phase diagram in the grand canoncial ensemble. 

5 Quark-antiquark potential 

We found above that a larger dilaton step-size, A, makes a baryonic phase more stable, 
which implies that a large A triggers a confinement transition as well as chiral symmetry 
breaking. Therefore, it is interesting to study the confmement/deconnnement transition 
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property by other methods. One of the standard tools is the Wilson-Polyakov loop, which 
can be obtained holographically by computing the on-shell Nambu-Goto action[3, 4, 39-41]. 



^ = i/drWl+f;) 4 /(^, (5.1) 



2ira' J V \RJ 

where the static string worldsheet is parametrized as x(r) and x is one of the field theory 
spacial direction. There are two types of solutions: a pair of parallel strings and a U-shaped 
string. In principle, the string end points can end on the deep IR bottom of the D7 probe 
brane. Such a non-trivial embedding is determined by the DBI action. However, in this 
paper, we will simply consider the D7 brane infinitely far away from the horizon, which 
corresponds to an infinitely heavy quark. 

For a pair of parallel strings, the embedding is simply 

x(r) = constant , (5-2) 

and for a U-shaped strings, we have the conserved quantity (x is a cyclic coordinate) 

y/&m + ijx* R R 

where ro is the minimum value of r giving x'{tq) = oo. From (5.3) we have 

x(r) = 2R 2 / dr = (5.4) 



" 2 \/7H^e^)r4//( e ^o)r 4 /(ro)) - 



We can compute the on-shell action of the U-shaped string embedding (S^ G ) by plug- 
ging (5.4) into (5.1). However, it is divergent because of the infinitely heavy quark mass 
contribution. We regularize the on-shell action by subtracting this infinite mass, which is 
nothing but the on-shell action of a pair of parallel strings (S^ G ). 
In summary, we define the quark-antiquark potential (V qg -(l)) as 

1/ n\ = ^ng ~ ^ng 
Vq " [L) - 2vra' 

f°° e'KO/ 2 f°° ,. . .„ ( 5 - 5 ) 

= 2 / rW 2 / dre^)/ 2 . 

Jro Jl-e*Wr 4 /(r )/(e*(r)r 4 /(r)) Jr a 

which is the function of I = x(oo), the distance between the quark and antiquark. ro on 
the right hand side can be related to I = cc (oo) by (5.4). Our numerical plots of V qq - for 
given parameter sets are shown in Fig 7. 

Let us start with the case A = which is just the well known J\f = 4 theory, Fig 7(a). 
This is the case with a trivial dilaton, = 1, and studied first in [3, 4, 39, 40]. At zero 
T, the potential scales as l/l (dotted line), which is a consequence of the conformality. At 
finite low T, the potential scales as 1/Z(1 + C((T/) 4 )). Again the dependence on the Tl 
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(a) A — 0, wh = (dotted line) and wh = (b) wh = 0.5, A = 0,1, 3, 5, 10 from top 
0.5 (solid line) 




(c) A = 3, w H = 2, 1, 0.5, 0.2 from top (d) A = 10, w H = 2, 1, 0.5, 0.2 from top 



Figure 7. Quark-antiquark potential (r = 1,A = 1.715) 

combination is due to the underlying conformal symmetry. The new feature at finite T 
is the existence of the phase transition from the bound quark antiquark pair to the free 
quark state, as the distance between quarks increases. In Fig 7(a) it corresponds roughly 
to the transition from the blue solid curve (a U-shape string) to the red horizontal line (a 
pair of parallel strings) near I = 1.5 2 . 

Now let us turn on A at a fixed T, Fig 7(b) {wh = 0.5 and here again we set the 
scale in (2.4) A = 1.715). As A increases a linear potential starts forming, which is the 
characteristic feature of confinement. It is interesting that the linearity is very clear at 
A = 10, where there exists a stable baryon phase Fig 5(b) 3 . So, we see a rough correlation 
between our phase diagrams Fig 5 and V q q (Fig 7). Fig 7(c) corresponds to Fig 5(a). The 
linearity is not clear in Fig 7(c) and it is natural that we expect a deconfined phase and no 
stable baryons as shown in Fig 5(a). Fig 7(d) corresponds to Fig 5(b). At high temperature 
we tend to loose the linearity and there is no big difference from Fig 7(c) (deconfinment, 
no baryon). However at lower temperature, the linearity becomes stronger (confinement, 
baryon phase). So it is consistent with the phase diagram Fig 5(b). 

2 In principle one should also take into account "graviton exchange" between the separated string world- 
sheets as was done carefully in [42] for large quark separation. This will modify the potential at large 
separation and will replace the phase transition shown at I ~ 1.5 with a cross-over. Since we focus on the 
dilaton effect at intermediate distance we don't consider that modification. If we considered it we could 
have a smoother cross over at larger distance, and also possibly, lose the additional first order transition in 
the inset of Fig 7(d) . 

3 This linear potential was observed also in a similar setup with — 1 — q/r% log (l — rjj/r 4 ^, where q 
plays a role of our A [41]. 



There is one additional interesting feature at low temperature and large A. As shown 
in the insert of Fig 7(d), there is a first order phase transition between the linear potential 
and Coulomb-like potential. This is at odds with QCD, but makes good sense in our model. 
Our dilaton profile is a step function so there is a potential barrier localized at r = A. As 
I increases, the string will extend inside, reduce ro and finally meet ro = A, where a linear 
potential is built. If we further increase Z, then the string finally will manage to get through 
the barrier at r = A and beyond that point, there is no friction against the string moving to 
lower r since the dilaton is constant. Therefore the potential becomes Coulomb-like again. 
We don't see this first order transition at high T, because at high T, the interval between 
th and A is too short and there is no room for a "Coulomb" phase, i.e. After getting inside 
r = A, the string should meet the black hole horizon very quickly. 

The Coulomb-phase should not be present in QCD. We could easily cure this artifact 
by considering a continuously growing dilaton profile in the IR rather than the flat form. 
However, to keep the phase diagram seen in Fig 5(b), the dilaton should only increase 
mildly (For example, as l/r q as studied in [43]). If it is too strong, the phase diagram will 
be always the type of Fig 5(c), destroying our motivation for this paper. However, a simpler 
explanation of the Coulomb phase is that it is a result of the quenched approximation. In 
reality we would expect qq pair production to break the string at separations before the 
Coulomb phase sets in for large A. For this reason these models still seem reasonable 
for QCD. In this section, our discussion has only been qualitative. More quantitative 
studies of the Wilson-Polyakov loop, together with other modifications of dilaton, would 
be interesting. 

6 Summary 

In this paper we have used the D3/D7 system to holographically study the phase diagram 
of a chiral symmetry breaking gauge theory as a function of the running coupling profile. 
We have included a running coupling through a phenomenological non-backreacted dilaton 
profile which steps between conformal UV and IR regimes. Here we have considered depen- 
dence on the height of the step. The model has previously been shown to have three phases 
[13]: a chirally symmetric quark plasma at high T and fx; a chiral symmetry broken phase 
at small T and /i; and a more exotic chiral symmetry broken phase with quark density at 
intermediate T, [i for some parameter values. The order of these phase transitions depend 
on the height of the step. A first order transition to the chirally symmetric phase can be 
achieved for low step values. Previous work has also shown that the low n transition with 
T can be made second order by phenomenologically tinkering with the shape of the AdS 
black hole horizon (in a way compatible with the spatial symmetry breaking of the D3/D7 
system). 

The crucial extra ingredient we have concentrated in this analysis is the low T and 
intermediate [i baryonic phase. Baryons can be introduced as D5 branes wrapped on the 
5 5 of the dual geometry and linked D7-D5 systems describe the hadronic phase. In [14] 
such configurations were introduced in a geometry with a diverging dilaton in the IR. The 
resulting baryonic phase persisted though to arbitrarily high density unlike in QCD (but 
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in a similar fashion to the equivalent phase diagram in the Sakai Sugimoto model [37]). 
Our intuition for our analysis was that in a model with a step function dilaton profile 
there would be no very large or small baryon vertices because in the two conformal regimes 
they would shrink away. The only D5 brane configurations would lie around the step 
and we could hope they would only play a role at intermediate \i. We indeed find, after 
careful numerical analysis, that this is the case for some range of the step size. Our model 
which best matches expectations for QCD in the low T regime is shown in Fig 5(b). The 
second order transition to this regime matches expectations in the absence of internuclear 
interactions which we neglected. 

We have also made a qualitative link between the presence of the baryonic phase and 
the observation of linear confinement in Wilson loop computations in the background. The 
step function dilaton form generates both confinement and chiral symmetry breaking when 
the step size is large enough. 

The phase diagrams we produce are not intended as predictions for QCD since the 
underlying physics model is somewhat different (eg the presence of super-partners) but they 
do demonstrate the wealth of behaviours possible in strongly coupled gauge theory. We 
can also hope for some universality and by finding models that match QCD's expectations 
for the phase structure it may be possible in the future to study phenomena beyond that 
stucture using the models. For example in [27] the temporal behaviour of these systems, 
such as bubble formation, could be followed through a first order phase transition in a 
model of this type. 
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